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Abstract 

We report results from a study of the charmed double strange baryons $7^ and Q*^ 
at Belle. The is reconstructed using the Q~tt'^ decay mode, and its mass 

is measured to be (2693.6 ±0.3+};5) MeV/c^. The nf baryon is reconstructed in 

the mode. The mass difference M^.o — M^o is measured to be ^70. 7 ± 0.9^Qg^ 

MeV/c^. The analysis is performed using 673 fb~^ of data on and near the T{4S) 
collected with the Belle detector at the KEKB asymmetric-energy e^e~ collider. 



1 Introduction 



The experimental study of charmed baryons, including measurement of their 
masses, widths and decay modes, is an important test of many theoretical mod- 
els that provide quantitative predictions for the properties of heavy hadrons. 

Unlike the A^, and even their excited states, a thorough exper- 

imental study of the charmed double strange baryon is long overdue. The 
("^^ ~ (i)^) is heaviest known singly charmed hadron that decays 
weakly. The quark content of the is c{ss}, where the ss pair is in a sym- 
metric state. There are many theoretical models that predict the mass of 
the However, the range of these predictions is rather wide: 2610-2786 
MeV/c^ pp. Among several published measurements of the mass only 
three report a statistically significant signal: E687 obtains (2699.9 ± 1.5 ± 2.5) 
MeV/c^ in the T,~^K~K'~tt'^ channel [2], using five decay modes CLEO 

reports (2694.6 ± 2.6 ± 1.9) MeV/c^ [3], BaBar quotes (2693.3 ± 0.6) MeV/c^ 
in their analysis of Q*^ — > f2°7 [1]. The world-average value of the mass is 
(2697.5 ±2.6) MeV/c^ 0. 
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The decay mode f2~7r+ was observed for the first time by E687 [6]. 

They reported 10.3 ± 3.9 events over the background of 5.8 events. With 13.7 
fb~^ of data CLEO also found a signal of 13.3 ±4.1 events in this channel [3]. 
FOCUS observes 23 ± 7 events of this decay mode [7]. With a data sample of 
230.7 fb"^ BaBar reports a raw signal of 156 ± 15 events [1]. 

The spin excited = (|)+ fif state was observed by BaBar only recently 
They measured the difference between the masses of the VL*^ and to be 
(70.8 ± 1.0 ± 1.1) MeV/c^. Different QCD-based models yield a rather wide 
range of predictions for the mass difference, from 50 MeV/c^ to 70 MeV/c^ [5]. 
A lattice QCD calculation finds the value (94 ± 10) MeV/c^ 0. 

In this work we present our results on the measurement of the fi^ mass using 
the decay fi^ VL~'k'^ [TD]. We also confirm the recent BaBar observation of 



2 Data and selection criteria 

The data used for this analysis were taken on the T(4S') resonance and in the 
nearby continuum using the Belle detector at the e^e~ asymmetric-energy 
collider KEKB [11]. The integrated luminosity of the sample is equal to 605 
fb~^ of on-resonance data and 68 fb^^ of data below the T(45'). 

Belle is a general-purpose detector based on a 1.5 T superconducting solenoid; 
a detailed description can be found elsewhere [12]. Tracking is performed by 
a silicon vertex detector (SVD) composed of concentric layers of double-sided 
silicon strip detectors, and a 50 layer drift chamber (CDC). Particle identi- 
fication for charged hadrons, important for the measurement of final states 
with kaons and/or protons, is based on the combination of energy loss mea- 
surements (dE/dx) in the drift chamber, time-of-fiight measurements (TOF) 
and aerogel Cherenkov counter (ACC) information. For each charged particle, 
measurements from these three subdetectors are combined to form likelihood 
ratios in the range from to 1, 

P{A/B) = CiA)/{CiA) + C{B)), 

where C{A) and C{B) are the likelihood values assigned to the K,7r or p 
identification hypothesis for a given track. From now on we use these likelihood 
ratios to identify protons, charged kaons and pions with criteria that have 
efficiencies greater than 85% for p/K/n. The probability for each particle 
species to be misidentified as one of the other two is less than 12%. 

An electromagnetic calorimeter (ECL) is used to detect photons, which are 
reconstructed from isolated ECL clusters that have no corresponding charged 
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M(AK")-M(A)+m^ [GeV/c 



Fig. 1. Q- ^ AK': M{AK-) - M{A) + tua spectrum of the selected Ag- 
nations. The selection requirements and the fit are described in the text. 



combi- 



track, and a shower shape that is consistent with that of a photon. 

A hyperons are reconstructed in the A p7r~ decay mode, fitting the p and 
vr tracks to a common vertex and requiring an invariant mass in a ±4 MeV / 

4cr) interval around the nominal A mass value [5]. We then apply the 
following cuts on the A decay vertex: 

• the difference Aza in the z-coordinate between the proton and pion at the 
decay vertex should satisfy A^a <2 cm, where the z-axis is parallel to the 
e"*" beam direction; 

• the distance between the vertex position and interaction point (IP) in the 
plane transverse to the beam direction (Ar^) should be greater than 0.5 
cm; 

• the angle l3\, between the A momentum vector and the vector pointing from 
the IP to the A decay vertex, should satisfy cos Pa > 0.0; 

• the vertex fit should have an acceptable x^- 

We apply only a weak cut on /3a since the analysis relies on A's emerging from 
a displaced vertex rather than from the IP. 

A sample of hyperons is then reconstructed in the decay mode 
AK~ . The K~ does not originate from the IP; we therefore impose only a 
loose requirement on the impact parameter of the kaon track, requiring it to 
be larger than 0.01 mm. We then require the difference between the AK~ 
invariant mass and the A mass to be in the range 1669 MeV/c^ < M{AK~) — 
M(A) + ruA < 1677 MeV/c^ (±3cr), where rriA is the nominal A mass [5], fit 
the A and the K~ track to a common vertex, and apply the following cuts: 
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• the distance in the transverse plane between the interaction point and the 
vertex position (Ar^^), obtained by fitting the A trajectory and the K~ 
track to a common point, should be greater than 0.01 cm; 

• the distance from the IP to the decay vertex for the Vt~ should be less than 
the corresponding distance for the A; 

• the VL~ momentum vector should point to the IP, satisfying cos a^^- > 0.99, 
where the angle a^- is measured between the Q~ momentum vector and 
the vector pointing from the Q~ decay vertex to the IP; 

• the angle a\, between the A momentum vector and the vector pointing from 
the Q~ to the A decay vertex, should satisfy cos > 0.99; 

• the Q~ vertex fit should have an acceptable x^- 

Figure [1] shows the M{AK~) — M(A) +mA distribution for the Q~ candidates 
after the requirements described above. A fit to this distribution with a double 
Gaussian to describe the signal and a third-order polynomial to describe the 
background gives 33880 ± 1680 Q~ AK~ events. The fitted Q~ mass of 
(1672.363 ± 0.012 (stat.)) MeV/c^ is in good agreement with the PDG value 
0. 



3 17° mass measurement 

Using the Q~ sample we reconstruct baryons in the decay mode — > 
without any requirement on their momenta to obtain candidates from 
fragmentation of cc quarks as well as from 5-meson decays. 

Figure [2|^a) shows the M{Q~'k~^) — M{Q~) + m^- distribution for the 
candidates: a clear signal peak is seen near 2700 MeV/ c^. The dotted histogram 
in Fig. [2]^a) shows the contribution from the Q~ sidebands, which is featureless. 
A fit {x^ /n.d.f. = 57.20/61) to the distribution with a Gaussian for the signal 
and a first-order polynomial for the background yields 725 ± 45 events at a 
mass of (2693.6 ± 0.3) MeV/c^; the Gaussian width was found to be (4.9 ± 0.3) 
MeV/c^. We exclude the region to the left of 2580 MeV/c^ from the fit to 
avoid a contribution from — > Q~7r~^7i^, where the tt^ is not reconstructed. 
To check the signal, we also reconstruct the wrong-sign combination Q^n^ 
(see Fig. [2](b)), where no peaking structures are observed, as expected. 

The systematic uncertainty on the measurement is assessed as follows. First we 
fit the signal in bins of Ar^^, Ar^, cosa^-, cosoa, p*{^c) p(^^)) where 
p*{Q^) is the reconstructed momentum of the candidate in the e~^e~ center- 
of-mass frame and p(vr"'") is the vr"*" reconstructed momentum in the laboratory 
frame. We divide the signal events into three roughly equally populated bins 
of the selected variables. In each pair of variables (Ar^^ and Ar^^, cosa^- 
and cosq!a, P*(^c) p(7r^)) we take the maximum value of the change of 
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Fig. 2. (a): n'^ 0~7r+: M{n^7r+) - M{n-) +mQ- spectrum of the selected $7"7r+ 
combinations. The selection requirements and the fit are described in the text. The 
dotted histogram shows the 0~ sidebands, (b): M{Q~TT~) — M{Q~)+'mQ~^ spectrum 
for wrong sign 0,^tt~ combinations. 

the fitted mass summed in quadrature, finding ^^'^ MeV/c^. We then vary the 
order of the polynomial in the fit function, fitting ranges and the Gaussian 
width within its error, finding no significant change in the fitted mass. Finally, 
since we actually measure M(f2°) — M{Q~) + m^-, the contribution from 
the f2~ mass uncertainty is ±0.29 MeV/c^ 0. All the sources of systematic 
uncertainty are summarized in Table [H We obtain: 



no 



2693.6 ±0.3 (stat.) (syst.)) MeV/c 



Our measurement of the ^2° mass is in agreement with the world-average 
value [5]. 



4 Confirmation of tiie f2*° 



From the sample we reconstruct the first excited baryon in the decay 
mode nf fi°7. We apply a 2679 MeV/c^ < M(fi°) < 2709 MeV/c^ re- 
quirement on the fl^ candidate. To suppress the combinatorial background, 
we require the photon candidate to have energy > 100 MeV and Qf can- 
didates to have momenta p* > 2.5 GeV/c, where is the energy detected in 
the electromagnetic calorimeter and p* is the reconstructed momentum of the 
Q*^ candidate in the e~^e~ center of mass. 

Figure [3](a) shows the mass difference M((]°7) - M(fi°) distribut ion for the 
remaining ^2^7 candidates: a clear signal peak in the mass difference M(f2°7) — 
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M{^,y)-M{Q.,) [GeV/c ] 




0.1 „ 0.15 „ 0.2 , 
M(D"y)-M(D") [GeW/c] 



Fig. 3. (a): Qf J]°7 (0° ^ J]-7r+): M(O07) - M(nO) spectrum of the selected 
fl^7 combinations. The dotted histogram shows the sidebands, (b): D*^ — > D^^ 
— > i^~7r"^): M{D^^) — M{D^) spectrum of the selected D^^y combinations. The 
selection requirements for (a) and (b) as well as the fit are described in the text. 

M{Q^) is seen near 70 MeV/c^. The dotted histogram in Fig. [3](a) represents 
the fl^ sidebands, which are featureless. A fit {x^/n.d.f. = 70.69/73) to the 
distribution with a Crystal Ball function [13] with parameters fixed from the 
Monte Carlo simulation for the signal and a second-order polynomial with 
an arctangent threshold term for the background yields 54 ± 9 events at a 
mass difference of (70.7 ± 0.9) MeV/c^. The fit gives a significance of GAa 
for the signal including systematics from varying the signal and background 
par ameterizat ions . 

The systematic uncertainty on the M{fl^'y) — M{Q^) measurement is assessed 
as follows. We vary the order of the polynomial and the threshold term in 
the fit function and signal width within its error, finding a mass shift of +0.1 
MeV/c2 and MeV/c^ respectively (Tabled. 



We study a possible bias in the measurement from the incomplete calibration 
of detector responses (energy scale and line shape) for low energy photons. 
Comparing Monte Carlo generated and reconstructed masses for Q*^ we find 
a shift in the mass of +0.9 MeV/c^. As a cross check, we reconstruct D*^ 
D°7, — > K~7r~^ using the full data set. Besides the particle identification 
selection described in Section 2, we require the D^-meson invariant mass to be 
within ±10 MeV/c^ of the nominal value (~ ±2(t), the photon to have energy 
> 100 MeV, and the D*^ candidate to have momentum p* > 2.5 GeV/c. 
By analysing the spectrum from the 50 MeV/c^ < M(^]°7) - M(fi°) < 90 
MeV/c^ mass window we find 82% of all 7's are in the energy range from 
100 to 200 MeV. Therefore we require the photons from the D*^ decays to 
lie in the same energy range. Figure [Hb shows the M(D°7) — M{D^) mass 
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Table 1: Contributions to the systematic uncertainty. 





Moo, 


Mn*o-M^o, 




MeV/c2 


MeV/c2 


Fit in bins 


+1.8 
-1.5 


- 




±0.3 


- 


Calibration mode 




-0.9 


Signal width 


0.0 


+0.1 
-0.2 


Fit function 


0.0 


+0.1 


Total 


+1.8 
-1.5 


+0.1 
-0.9 



distribution for the remaining D°7 candidates after these requirements. From 
a fit by a Crystal Ball function for the signal and a third-order polynomial for 
the background we obtain M{D^-f) - M{D^) = (143.03 ± 0.09) MeV/c^. This 
is 0.91 MeV/c^ higher than the world-average mass of m£)*o — rriDO [5], which 
is consistent with the Monte Carlo study. We assign a conservative systematic 
error of —0.9 MeV/c^ due to the uncertainty of calibrations. 

Combining these results (see Table [1]), we find 

Mn*o-Mno = (70.7 ±0.9 (stat.) ^ol (syst.)) MeV/cl 
in good ag reement with (70.8 ± 1.0 ± 1.1) MeV/c^ obtained in Ref. [1]. 



5 Summary 

In conclusion, we report a precise measurement of the mass using its decay 
to fi-7r+, Mno = (2693.6 ± 0.3 (stat.) +J;^ (syst.)) MeV/c^. The presented 
analysis provides a measurement of the 1)° mass based on the largest avail- 
able statistics. The result agrees with the world-average value [3] and has a 
significantly smaller uncertainty. 

We also search for the Q'f meson decaying into fic7 that was recently observed 
by BaBar [1]. We use our sample of fi^ — ^ Q~7r~^ decays and confirm the pro- 
duction of in cc fragmentation with a significance of 6Aa. The mass differ- 
ence M^,o - M^o is measured to be (70.7 ± 0.9 (stat.) (syst.)) MeV/c^ 
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